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ABSTRACT 

We re-analyze the M31 microlensing event WeCAPP-GLl/Point-AGAPE-S3 taking into account 
that stars are not point-Uke but extended. We show that the finite size of stars can dramatically 
change the self-lensing event rate and (less dramatically) also the halo lensing event rate, if events 
are as bright as WeCAPP-GLl. The brightness of the brightest events mostly depends on the source 
sizes and fluxes and on the distance distribution of sources and lenses and therefore can be used as 
a sensitive discriminator between halo-lensing and self-lensing events, provided the stellar population 
mix of source stars is known well enough. Using a realistic model for the 3D-light distribution, stellar 
population and extinction of M31, we show that an event like WeCAPP-GLl is very unlikely to be 
caused by self-lensing. In the entire WeCAPP-field (17.2' x 17.2' centered on the bulge) we expect only 
one self-lensing event every 49 years with the approximate parameters of WeCAPP-GLl (full-width- 
half-maximum (FWHM) time-scale between 1 and 3 days and a flux excess of 19.0 mag or larger in 
R). On the other hand, if we assume only 20% of the dark halos of M31 and the Milky- Way consist of 
1 solar mass MACHOs (Massive Astrophysical Compact Halo Objects) an event like WeCAPP-GLl 
would occur every 10 years. Further more, if one uses position, FWHM time scale, flux excess and 
color of WeCAPP-GLl, self-lensing is even 13 times less likely than lensing by a MACHO, if MACHOs 
contribute 20% to the total halo mass and have masses in the range of 0.1 to 4 solar masses. We also 
demonstrate that (i) the brightness distribution of events in general is a good discriminator between 
self and halo lensing (ii) the time-scale distribution is a good discriminator if the MACHO mass is 
larger than 1 solar masses. Future surveys of M31 like PAndromeda (Pan-STARRS 1) should be able 
to provide many more such events within the next 4 years. 

Subject headings: dark matter — gravitational lensing — galaxies: halos — galaxies: individual (M31, 
NGC 224) — Galaxy: halo — galaxies: luminosity function, mass function 



1. INTRODUCTION 

Microlensing searches towards Local Group galaxies arc interesting in at least two respects. They can constrain the 
fraction of compact halo dark matter (MACHOs, see Paczyiiski (1986); Gricst (1991)) and they allow to study stellar 
populations and the 3D distribution of stars in the Milky Way and the target galaxies. Lensing events caused by stars 
(self-lensing) also define a lower limit to the number of lensing events that have to be identified in a survey and therefore 
provide a consistency check for the lens model (depending on stellar popul a tion c ontent, stellar dynar nics and density 
distribution of the stars) and for the survey efficiency (see 'Alcoc k a/] (pOOT) for MACHO, and [T isserand et al\ 
([200l : lAfo^so et a l. (2003) for EROS and Calchi Novati et al. (20Q1) for POINT- AGAPE and Ide Jon g et al. (200^ 
for MEGA). One can use the known characteristics of lensing and self lensing events to design surveys that will be 
do minated by self lensin g. These self lensing-surveys can measure the faint end mass function of stellar populations 
(see iRiffeser et al\ (l2006h ). 

One can obtain the most likely MACHO-mass fraction and its confidence limits from the analysis of all lensing events 
found in a survey (using number, spatial distribution, time scale distribution etc.) only after the selection criteria 
and the survey efficiency have been taken into account. The survey efficiency depends on the event characteristics 
(location, color, time scale and flux excess, finite source effects etc.) and on the sampling and photometric quality of 
observations. Including precise values for the survey efficiency can completely change the interpretation of a survey. 
iPaulin-Henriksson fc Calchi Novatil (|2004l ) concluded that there is no ev idence for MACHOs towards M31 in the INT 
data set, whereas the same collaboration claimed with lCalchi Novati et a l. (2005) that there is a fairly strong evidence 
based on new efficiency estimates of the same survey. 

Instead of comparing the expected and observed events one can analyze the observables of individual lensing events; 
these are flux excess at maximum magnification, full width hal f maxi mum time, color, location and presence/absence 
of finite source signatures in the light curve (see iRiffeser et all (120061 ) for more details). One can ask for the relative 
probability of halo lensing and self lensing and derive pro bability distrib utions for the lens masses causing that ev ent. 
This has been done in a simphfied way for WeCAPP-GLl (jRiffeser et aU (2003), or Paulin-Henriks son etaP (l2003l) for 
the P0INT-AGAP E-S3 identification of the same event) and also for POINT- AGAPE Nl (Auriere et aLN200H) and 
N2 (jAn et aLll2004f ) in the past. However, the sources for all three events have been treated as point sources. This 
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Fig. 1.— Light curve of WeCAPP GLl (R band: blue, I band transformed to R band: red); POINT- AGAPE-S3 (r' band: dark green, i' 
band: light green) 

can mislead the interpretation, in particular if the events are very bright. This is also the case for light curves not 
showing any finite source signatures because the number of possible lens-source configurations strongly changes for 
bright events and the interpretation of the lens mass is modified. 

It is the subject of our paper to (re)-analyze WeCAPP-GLl and derive the relative probabilities of self lensing and 
halo lensing accounting for the extended sizes of stars. We summarize our M31 model in Appendix ^1X1 and present the 
equations needed for extended source stars in Appendix [Bl and ICl We demonstrate the pronounced differences resulting 
from this more accurate description of M31 in Sj2l In f|3]we present the lens mass distribution without accounting 
for finite source effects and in 21 show some qualitative results using the radius for the source stars. We evaluate 
in ^J5] an accurate estimate for the mass distribution of the event using the stellar source size. We investigate in ^ 
how sensitive the results are with respect to the mean M31 extinction and the line-of-sight (LOS) extinction towards 
WeCAPP-GLl and in Sj7]with respect to the stellar population properties of the disk stars acting as sources. In fj8]we 
give an outlook to the statistical interpretation for microlensing events in M31. Finally we draw our conclusions in fJ5) 

2. AN IMPROVED WECAPP-GLl ANALYSIS IN THE POINT SOURCE APPROXIMATION 

For the current analysis we improve our light curve fit for WeCAPP-GLl (Fig. [T]) with respect to iRiffeser et oil 
((200l . Table 2 by using a Paczyhski microlensing light curve fit (5-|-2-f 2-1-2 free parameters in 4 filters; to, uq, ^e, 
Fo,R, C'ri ^o,Ij Ci) Fo,t', Cr' , ^0,1' I C'i') for the observables and a Gould microlensing function to determine the errors 
of the observables ipwHM and AF. We also improved our photometric calibration accounting for the color terms in 
the filter calibration; our instrumental magnitudes R and / [phot /sec at AM=0] transform to the Johnson system with 
R = R + 23.58 - 0.01(^ - /) and / = I + 22.89 + Q.22{R - I). The resuhing parameters are listed in Table 1. Note 
that the fit implies that the source star of GLl is a bright Red Giant Branch (RGB) star. 
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TABLE 1 

Summary of WeCAPP-GLI (or AGAPE-S3) observables (a and <5) and parameters derived from a Paczynski microlensing 

LIGHT CURVE FIT: DATE OF LIGHT CURVE MAXIMUM, FULL- WIDTH-HALF-MAXIMUM TIME tpWHIvI OF THE EVENT, FLUX EXCESS AT LIGHT 
CURVE MAXIMUM IN R AND I-BAND, COLOR OF EVENT, QUALITY OF THE FIT, AND AN ESTIMATE FOR THE DEGENERATE AMPLIFICATION AND 

ABSOLUTE BRIGHTNESS OF THE SOURCE WITHOUT ERRORS. 



Despite the well-sampled light curve of the event and the nicely fitting Paczynski light curve, finite source signatures 
above the highest data point of 8 x 10~^ Jy cannot be ruled out. IVIoreover, as we will demonstrate later on, finite source 
signatures for events brighter than 19 mag are indeed mor e likel y than no finite source signatures. In a future work 
we are planning to use a modified light curve (see Riffes er et al\ ([20 06). Eq. 12) that allows to include the finiteness 
of lenses and sources and their limb darkening. 

The mass probability function for WeCAPP-GLl was already derived in IRiffeser et al\ ()2003f) using a simplified 
analysis of the event and a simple description of M31. 

We now use the improved analysis method derived in IRiffeser et al\ (120061 ) and a more detailed 1VI31 model. We use 
this 1VI31 model for all calculations further on unless stated otherwise. The relevant differences for the analysis of 
WeCAPP-GLl are as follows: 
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Fig. 2. — Left panel: Color (R — I) - _R-band magnitude relations for stellar populations of the bulge (red) and disk components (blue) 
with solar metalicities (using isochrones ofjGirardi et al. (2002) in the Johnson filter system). Right panel: Luminosity function gfor the 
bulge (red) and disk (blue) combining isochrones, weights and mass functions of our model. The green curve shows a nearly continuous 
disk population with 71 components. 

1. C olor-magnitude rel ation of source stars 

In i Riffeser erall (HQQI) we approximated the color-magnitude-relation pcmds{-M,C) of bulge and disk stars in 
M31 with observations of M32 to derive a brightness estimate for a post main sequence star with a color of the 
WeCAPP-GLl event. In this way, the brightness of a bulge star and a disk star with a color of WeCAPP-GLl 
was estimated the same, Mn — —1.5 mag. 

We no w model the disk and bulge population separately, using the stellar population isochrones of lGirardi et al\ 
([2002) with a metalicity of Z = Zq (isoc_z019.dat). We describe the bulge as a 12.6 Gyr single stellar 
population (SSP). 

The disk is modeled as a composite of 6 SSPs with ages of 4, 20, 100, 500 Myrs and 2.5 and 12.6 Gyrs. Their 
relative weights are calculated by integrating an exponentially declining star formation rate with t — 8 Gyrs over 
adjacent intervals (see Table 2). This simple model is close to a population where stars formed continuously 
over the past 12 Gyrs. In Figs. [Hand [3] we compare the luminosity and mass function of this composite stellar 
population with a population composed of 71 bursts with the same star formation decline rate of 8 Gyrs. The 
differences are marginal, i.e., our 'simple' model describes the brightness and mass distribution of stars with an 
e-folding timescale of 8 Gyrs fairly well, and offers the advantage of much faster numerical integrations in the 
color-magnitude plane. 



age [Myr] 


interval [Myr] 


weight 


4.0 


0.0 - 12.0 


0.00039 


20.0 


12.0 - 60.0 


0.00158 


100.0 


60.0 - 300.6 


0.00804 


501.2 


300.6 - 1506.5 


0.04417 


2511.9 


1506.5 - 7550.6 


0.35630 


12589.3 


7550.6 - 12589.3 


0.58952 



TABLE 2 



The star formation history of M31 is not known exactly^. Our choice ensures the presence of mostly old stars 
with a small fraction of young, blue and bright stars. 

This implies that the stars in our model also populate the extremes of the M31-CMD (see Fig. [5]). 

We thus can allow for the effect that a few bright stars might be more efficient sources for detectable microlensing 
events than many fainter, older stars. 

We now do not have to estimate the typical stellar brightness at a given color anymore, but can construct a 
brightness-color probability distribution Pcmd(A^,C) from the isochrones and the assumed stellar mass function 
£,{M). The information provided in the Girardi isochrones also allows a straightforward extraction of the stellar 
radii R^,{A4,C), which is necessary for an accurate treatment of finite-source size effects. 

^ IBrown et al\ Il2003l )'s data and analysis suggest that (at least at a distance of lOkpc wh ere their data hav e been taken) about 
30% percent of star s (in mass) c ould be as young as 6 — 8 Gyr. Also, the metalicity is large fsee lMould &: K ristianI 119861) . [M/H] =-0.6 
IDurrell et all llT99l ). fFe/H]=-0.6 IRich et ali (11999) [Fe/H]=-0.4 ] and falls off only at more than 30 kpc IKalirai et aLI 12009) . 
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Note that the color information of the event helps to select a smaller set of sources and therefore the mass 
probability results from a smaller range of source brightness. Because of this additional information the mass 
probability distributions become narrower and more precise. 

A nalysis method 

In lRiffeser et al\ ()2003D we had estimated the approximate source brightness A4 from the color of the event C"^^^^^ 
and derived a magnification Aq to convert the observed full width half maximum time of the event tpwHM its 
Einstein-time scale « t^^"^M/VT2(Ao - 1) « <S?^"&m/VT2 Af;^"7{^^VcgalO-"-4[^(P™^''^'"°'°)+"4-43l}. 
Using descriptions for bulge/disk/halo densities p and velocities pvt, for the mass functions of stars £,{M) and 
a single mass MACHO halo we obtained the relative probabilities p{M) for self lensing and halo Icnsing (as a 
function of the MACHO mass) with 

OO -Dos g 

p{MM)- j Ps{D,s) J pi{Doi) p^„(^^^ (^^^ dDoidDos . (1) 



This means that the lens mass probability distribution p{M, Ie) was obtained using the (unobservable) maximum 
magnification Aq and Einstein-time instead of the observables flux-excess A™"***^ and full-width half maximum 
time scale ^p\^jj|^ • 

We have shown in lRiffeser et al\ ()2006f) (Eqs. 80 and 84) that the Einstein-time-impact-parameter degeneracy 
leads to coupled errors for these two quantities, and thus should be avoided. In the following we will always use 
the real observables fl ux excess Af and tim e-scale ipwHivi instead of maximum magnification and Einstein time. 
Using the notation of lRifi'eser et ali ()2006 '). the event rate per area, per event time scale tpwHivi, per flux excess 
AF, per source color C, per absolute magnitude of the source star M, and per lens mass M is 



d'^T{x, y, tpwHM, AF,C, M,M) ^ 2p^^d^ 
dx dy dtpwHM dAF dC dM dA4 tpwHM"^ 



J pR%Pvt{vt)dDoi + il* p* j Pv^{vl)T*'^ duQ 



dDos 
(2) 



All definitions are explained in Appendix [Bl 



3. Mass functi o n of stars 

iRiffeser et al\ (120031 ) assumed that both disk a nd bulge star s are c onfined to a mass interval from O.O8M0 
to 0.95Mf7) /V'^^se) and to lOAf© (disk) with a IZoccali et aU ()2000[ ) IMF (^ ~ M'^-^^) for the bulge and a 

iGould et al\ (|1997f ) IMF for the disk, respectively. W e now use tho s e mas s f unctions which are consistent with 
the st ellar populat ion models, i.e., a combination of a lBallero et al\ (|2007f ) fc iKrourial ([2002) IMF for the bulge 
and alGould et al] (1997,1 fc .Kroupal (|2002il IMF for the disk. Mass loss is provided by stellar population models 
(jGirardi et aLl 120021)" and is taken into account in the mass function (MF). 



We also extended the lower mass limits to O.OO17\f0. 

Because of the small gradient of the MF for very small masses and their low cross section for microlensing all 
results are independent from the lower mass limit. 

We al so included stellar remnants originating at the high mass end of the IMF into the MF fsee lRenzini fc Ciottil 

(EllD). 

4. Extincti on 

In Riffeser et al\ (|2003f ) we did not account for the extinction along the line of sight to the source of the lensing 
event. Extinction alters the true color and brightness and thus the necessary magnification of a lensed source, 
which changes the shape of the mass probability function extracted from individual events. 

The overall extinction decreases the observed flux, leading to the need to increase the number of stars (and thus 
the mass) to obtain the same luminosity. 

Extinctio n in this way a lso i nfluences the absolute and relative amplitudes of lensing and self lensing rates. The 
model of IRiffeser etal\ (|200al assumes an on average M31-extinction for all disk stars in the WeCAPP field of 
0.51 mag and for all bulge stars of 0.19 mag in the R-band, independent of the angular position of the event, and 
independent of the LOS distance to the source. The MW-extinction is added uniformly with 0.17 mag (R-band). 

With these assumptions on extinction one gets the following extincted colors for the bulge and disk stellar 
populations described above: for the bulge {U — i?)ineas — 0.78, {B — V^)mcas = 1-19, {B — i?)mcas — 1-91, 
(V - RUe.s = 0.72, (i? - /)„,eas = 0.75, and for the disk {U - B)„,eas = -0.06, {B - VU,,, = 0.73, {B - 
-R)meas = 1.33, (V — i? ) meas = 0.60, (i? — /)incas = 0.69. Especially the values for the disk agree well with 
IWalterbos fc Ke nnicuttI (|1987f ). We will use these extinction values if not specified otherwise. Additionally we 
use slightly refined extinction descriptions in § [6] later on. 
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Fig. 3. — Present day mass function (MF) for M31 bulge (red) and disk (blue) population (mass loss according to lGirardi et al] l|2002l)). 
Stars more massive than IMq are either young disk stars (blue) or disk/bulge stellar remnants (white dwarfs < I.ISMq, neutron stars 
?s IAMq, black holes > 20Mq). The green curve shows a nearly continuous disk population with 71 components. 

5. Mass-to-light ratio [M/L] and total mass of bulge, disk and halo 

Combining the luminosity function and the mass function results in the M/L ratios of {M/L)fi — 3.3 (bulge) 
and {M/L)ii — 1.2 (disk) not including extinction. 

Using these (M / L) with the M31-light-model and the extinction assumptions results in a very similar massive 
bulge (4.4 X 10 M0) and a more massive disk (4.2 x IQi^'^Mq). Since we changed in particular the disk mass 
we now use a different halo model with a core radius of 5 kpc, a total mass of 1.23 x lO^^M©, and a cut-off 
radius of 100 kpc. Note that the halo model has a large uncertainty since a large set of different combinations 
of core-radius, central-density, cut-off radius are able to fit the measured rotation curve of M31. 

3. POINT SOURCE APPROXIMATION 

In this section we present the lens mass distribution without accounting for finite source effects. We no w demonstrate 
the im pact of the change of assumptions described above on the interpretation of the event relative to iRiffeser et al\ 
(|2003f ). In this chapter the sources are still treated as a point source. 
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Fig. 4. — Lens mass distribution [event rate r(Af) per area, per event brightness A_F, per event time scale 



per color (R — I)\ 



and their errors (Eg. IC10|I . In the left panel the measured event color (C™<=!is) was transformed to the un-reddened color (-R — /)o = 0.59 
for disk sources, and (R — I)q = 0.70 for bulge sources, both with a Gaussian error of 0.03. The right panel was calculated not using the 
color information of the event. The red, green, blue and purple curves show the event rate for bulge-bulge, bulge-disk (bulge lens and disk 
source), disk-bulge and disk-disk lensing. The mass range is confined to the MF interval, and therefore there is an upper limit of about 
50Mq resulting from black hole remnants. The amplitudes of the lens mass distributions are scaled such, that their maxima represent 
the expected rate due to the different lensing and self lensing scenarios and the ratios of their maxima yield the relative probabilities for 
the different scenarios. The information about halo lensing (M31-halo-bulge and M31-halo-disk data points are in black and grey, and 
MW-halo-bulge and MW-halo-disk in brown and light-brown) is displayed as points and not by a curve: we assume a mono-mass-spectrum 
for the halos, and each point represents a halo which is made of a given MACHO-mass to 100%. For each MACHO mass one obtains the 
event rate from the amplitude at that mass. 
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The left panel of Figure [4] shows the lens mass distribution (event rate T{M) per area, per event brightness AF, per 
event time scale tpwHivi) for WeCAPP-GLl in the point source approximation using the LOS position (j^meas^ymeas^^ 
the time-scale (ipwHivi)' the flux-excess {Ap'^'^^) and the col or (C'^^'^) of the eve nt and their errors. The assumptions 
in evaluating the event rates have been changed relative to lRiffeser et al\ (j2003D as detailed in (1) to (5). The event 
rate for bulge-bulge lensing is 2 times as likely as bulge-disk lensing, about 13 times as likely as disk-bulge lensing and 
about 12 times as disk-disk lensing. For lensing of a bulge star by a 0.7AMq MACHO in M31 we obtained an event 
rate which is 12 times as likely as bulge-bulge self lensing. Using the event color the most likely halo lens masses for 
WeCAPP-GLl are between O.ISMq an d 5.6Mf7^. The main difference of this new analysis (still in the point source 
approximation!) relative to the result of lRiffeser et al\ (|2003f ) is tha t the M31-halo - disk l ensing scenario becomes now 
up to a factor 1/4 similar to the halo-bulge lensing (compared to iRiffeser et al\ (|2003'), where the M31-halo-bulge 
lensing was about 10 times as likely as halo-disk lensing). For the point source approximation the total halo-lensing 
contribution (summing over M31 and MW halo-bulge and halo-disk lensing) is 13 times more likely than all self lensing 
contributions. 

In the right panel of Figured we show the same calculations without using the color information of the event. At 
low masses the curves for T{M) are only changed for disk sources, whereas bulge sources are not affected. The high 
mass end cut-off however is suppressed for all lens-source configurations due to the use of the measured event color. 
This is because the color of the WeCAPP-GLl makes the source most efficient for becoming a very luminous event, 
and sources with colors different from WeCAPP-GLl have smaller fluxes requiring larger Einstein radii and thus larger 
lens masses to achieve an event of the same brightness. Therefore, dropping the color information for WeCAPP-GLl 
mainly affects the high mass probabilities. 

In Fig. Owe show the distribution of events in the Ai^-tpwHivi-pla-ne in the point source approximation as obtained 
from Eq. [C8]with location and intrinsic color of WeCAPP-GLl, [{R - I)o = 0.59 for disk sources, {R - I)a = 0.70 
for bulge sources, without using its error]. Since in the point source approximation arbitrary large magnifications are 
possible, there is no limit for the brightness or shortness of the events. The event rate for events like WeCAPP-GLl is 
equal for bulge-bulge and bulge-disk lensing. MACHO-lensing events caused by a 0.2 Mq MACHO halo are roughly a 
factor of 10 more likely than self lensing for a full MACHO-halo. 

Here the color information affected the shape of the event-rate contours in the tpwHM-Ai^-plane. The 'tilted M'- 
shape comes from the bimodal luminosity distribution of stars with color of WeCAPP-GLl: The l uminosity function 
maps the magnification-tFWHM-distribution into the AF space, see Eq. 62 in IRiffeser et al\ (j2006( l. The color of the 
event constrains the luminosity of the post main sequence (PMS) sources fairly well, and one thus gets the 'tilted 
M'-structure of the event rate in the AF-tpwuM-pla-ne. The left part is due to main sequence source stars, the right 
part due to PMS stars. Main sequence stars need larger magnification for the same event brightness, and since large 
magnification needs small impact parameters b, this goes in parallel with short FWHM-timescales, ^fwhm c)c b/vt, 
of the events. If the PMS sources have a broader brightness distribution, the right part of the 'M' is washed out by 
moving it up and down vertically. Integrating over all colors makes the 'M-shape' almost disappear, since the width 
of the 'M' depends on the color (sec Riffescr et al. (2006) for more details). Since the intrinsic color of WeCAPP-GLl 
for bulge sources is between 0.6 and 0.8 (depending on the degree of reddening, see chapter [6]), and the magnitude of 
the PMS-stars is fairly constant in this color range (see Fig. [5]), the location of the right part of the 'M' is independent 
from color within the possible range of intrinsic colors of WeCAPP-GLl. The location of the event relative to the 
'tilted M' provides an indication weather the time scale of the event fits to its brightness. 

Figure [5] can also be used to visualize how plausible the observables AF-^fwhm of an event candidate are. E.g., if 
MACHOs would have 100 solar masses, then for a given brightness of GLl, its timescale would be too low - or, on the 
other hand assuming GLl is a convincing event, for MACHOs with IOOMq one would expect many more events with 
the same brightness but timescales of 10 days (which are not found). 

To summarize in the point source approximation an event as bright as WeCAPP-GLl is quite unlikely to be a 
self lensing event. The most likely self lensing scenario is bulge-bulge self lensing. The color information does only 
marginally influence the mass distributions. Noticeable is that for self-lensing low lens masses of O.lAf0 achieve relative 
high probabilities. Self-lensing is half as likely as halo-lensing, if MACHOs contribute 20% to the total halo mass and 
have a mass of about IMq. 

4. EFFECTS OF FINITE SOURCE SIZES: A QUALITATIVE UNDERSTANDING OF WECAPP-GLl 

In this chapter we show for the fittings parameters of the source star (color {R — /)meas = 0.83 and brightness 
■A^moas — —0.64) that the lensing interpretation is strongly influenced by the finite-source treatment. Here we discuss 
the different competing effects qualitatively, while the full quantitative formalism is described in SjSl 

If a source is treated as point-like, it can be magnified at light curve maximum by a foreground point mass with any 
magnification value between one and infinity: for each Einstein-radius of a potential lens, there exists a source-lens- 
trajectory (an impact parameter) such that the desired magnification of the event can be achieved. However, if the 
source is extended, then the magnification at light curve maximum is no longer unlimited: 

• if the impact parameter is larger than the projected source radius [6 > 6* « i?*£'oi/(2I?os)], nearly no finite 
source signatures appear in the event light curve, and all relations for the point source approximation hold. 

• if the impact parameter is closer {b < b*), the maximum magnifica tion does only weakly depend on the impact 
parameter anymore. In our approximation fsee lRiffeser et all ()2006[ ). Fig. 1) the maximum magnification depends 
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Fig. 5. — Contour plot of the time-scale and flux-excess distribution (AF-tpwHM-plane) of events at the location and with the intri nsic 
color of WeCAPP-GLl, [{R — I)o = 0.59 for disk sources, (R — I)o = 0.70 for bulge sources, without using its error], obtained from Eg. IC8I 
in the point source approximation. The first row shows self lensing, the second row M31 halo-lensing, the third row MW halo-lensing. 
The two further observables of WeCAPP-GLl (flux excess and event time s cale obtained from a point source lensing light curve fit) are 
marked in green. The estimate for the WeCAPP 6a detection limit (see IRiffeser et ol] II2006I) . Table 1) at the position of GLl is shown 
as white line. The contours differ by factors of lO"'^, and the dashed contour is that of 10~^ events per year, per square arc minute at 
the location of WeCAPP-GLl, per color (in magnitudes), per logarithmic timescale (time scale in days) and per flux-excess (converted to 
magnitudes) (see footnotes [3] and [2]l . 

solely on the Einstein-radius and on the size of the source star, projected onto the lens plane, i^^DoiZ-Dos (see 
below, Eq. ([3])). Therefore, for the highest magnifications there is no longer a trade-ofi^ between source- lens 
impact parameter b and Einstein-radius Re, but between projected source size i?*Doi/I3os and Einstein-radius 
Rb- This means that the larger the projected source size is, the larger the Einstein-radius has to be to achieve 
a high magnification. 

For self lensing, the Einstein-radii are limited by the maximum stellar masses (e.g. around IMq for the M31 bulge 
population) or the maximum remnant masses, the inefficient lensing geometries (fairly similar source Dos and lens 
distances Doi), and large projected source radii. Therefore, high magnification events are hard to achieve. These 
effects are much less severe for halo lensing events. One therefore expects for bright lensing events that halo lensing 
becomes more likely relative to self lensing if the finite stellar sizes are taken into account. This could make the self 
lensing hypothesis for an event as bright as WeCAPP-GLl very unlikely, although not completely un-feasible. 

How easily an event can occur depends on the number of source-lens pairs that can produce the event. The LOS 
density distributions of bulge and disk are confined to a few kilo parsecs (see Fig. [51 left) which implies that the 
source-lens distances along the LOS must not be larger than a few kilo parsec for self lensing; otherwise either the 
source or the lens density and thus their product is small. 

The relation between ma ximum magnifi c ation = (A_F + Fq)/Fq w AF/Fq, Einstein-radius i?E, and projected 
source radius R^Doi/Dos is (jRiffeser et al\ pOOl . Eq. (18)) 
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Fig. 6. — Left pa nel Density of bulge (red) and disk (blue) stars along the LOS to WeCAPP-GLl, using the M31-model described in 
IRiffeser et al\ I I2006I) . In the middle panel we have transferred the stellar brightness to a flux (in Jansky's) at 770kpc (distance of M31), 
and obtained the ratio of this flux to the stellar size (in solar units). This ratio is an indicator of the maximum brightness of lensing events 
in the extended source description (see Eq. [Sjl . Right panel Minimum distance that a source lens pair must have in order to cause an event 
as bright as WeCAPP-GLl (A™"^^*^ = 10.15 X 10~^ Jy) as a function of the source characteristic (flux and radius) and lens mass. The light 
red, magenta, green, cyan, blue, black lines represent lenses with masses O.1Mq,O.5Mq,1M0,5Mq,1OMq, and 50Mq, respectively. 




Da 



where the un-lensed source flux Fq — fluxR(A1 + .4)(10pc/Z?os)^ (sec footnote [3l) depends on the absolute source 
luminosity A4, the source distance -Dos and the extinction A towards the source. Note that our approximation Eq. [3] 
represents also the largest flux excess achievable by a light curve not showing finite source signatures. 
Eq. [3] can be inverted to 

^ol max -Dos(l+( Af,nax ^'^^'^ ^ 7^^ ■ (4) 

Eq. [4] provides the maximum distance a lens can have (to us) to allow a lensing event with a flux-excess Ai^max at 
maximum magnification once the finite source sizes are taken into account. For a given source distance and a given 
source and lens population, the upper limit of that distance is set by the largest Fo//?*-ratio for source stars and the 
largest lens mass possible. For a given (measured) flux excess of an event, Eq. (|4]) depends on the extinction along the 
line of sight: the smaller the extinction, the larger is A-Fmax! also, the radii i?* have to be taken into account using the 
de-reddened color of the event. If the color information is available, one has to use the stellar radii of the de-reddened 
source to obtain the largest Fq/R* - ratio for source stars. 

In Fig.[6}D (middle panel) we show the Fo/i?*-ratios for a bulge 12.6 Gyr SSP and 5 younger components of the disk 
composite SP, where the flux _Fo is given in Jansky (Jy), the radius i?* in units of the solar radius, and the source 
is placed at 770 kpc. The brightest bulge events, that can occur, have an un-reddened color of about (i? — /) w 0.7. 
For the ages above 500 Gyr the stars which can cause the brightest lensing events are either very blue, or have a 
color of about {R — I) « 0.7 (i.e. PMS stars present for the ages shown here). The very red luminous stars with 
(i? — /) > 1 are ineffective sources for bright lensing events since they are so large, implying a smaller magnification 
at light curve maximum. The numerous (but fainter) turn-off stars of an 12.6 Gyr old stellar population would need 
a reddening of E{R — /) « 0.45 mag to show the measured color of (i? — /) = 0.83 mag. Because this translates into 
an extinction of Ar w E{R - /)/(! - {Ai/Av)/{Ar/Av)) = 0.45/(1 - 0.482/0.748) 1.27 mag these stars can only 
produce {R — I) = 0.83 events which are 13 times fainter than the brightest bulge stars. 

The position of GLl (3'6" south, 2'40" west of the M31 center) supports the assumption that the source star belongs 
to the bulge population. Using the assumed extinction for the bulge of Ar = 0.36 mag the intrinsic color of GLl 
changes to (i? — /)o = 0.70, which is also the source color for the brightest bulge events (see the peak of the fiux-to- 
radius relation for a 12.6 Gyr old bulge star in Fig. [Bb). Such a bulge star has an un-extincted absolute brightness 
of Mr, = —0.91 mag and a radius of i?* = 30i?o. This is in good agreement to the (extinction corrected) R-band 
brightness of -1.00 mag obtained from the light curve fit (see Table 1 ). 

In Fig. [51: (right panel) we show the minimum lens-source distances (I?os — F)oi) as a function of brightness, radius 
and mass Fo/{R*/Rq) x ^M/Mq using the measured flux excess of WeCAPP-GLl A^°^" = 10 x 10.07-^ Jy. For 
bulge stars the largest FolO~°-^-^/(i?*/i?0) values are of order 4 x 10^'* Jy (middle panel, red curve). Now using this 
with the flux excess of GLl combined with a bulge extinction of Ar — 0.36 mag, a 1M0-Iens needs a minimum source 
distance of about 8.1 kpc (green line in the right panel of Fig. [5]), a separation where the product of source and lens 
densities becomes small (Fig. |B^). 

In Fig.[7|we combined the Figs.[B^ andlHJ; for the LOS of WcCAPP-GLl and show the minimum source-lens distances 
as a function of the source distance for different self-lensing masses (Mq — 0.1, 0.5, 1, 5, 10, 50 M0). The two peaks 
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Fig. 7. — Minimum lens-source distance for all self lensing configurations for the WeCAPP-GLl event, using the intrinsic color of 
WeCAPP-GLl of {R — /)o = 0.70 for bulge sources and (R — I)o = 0.59 for disk sources and Eq. (|4)l. The red, magenta, green, cyan, 
blue, black lines represent lenses with masses O.1Mq,O.5Mq,1Mq,5M0,1OMq, and SOA/q, respectively. The contours show the product 
of the density of all source and lens stars (irrespective of color) as a function of the LOS distance to GLl and the lens-source distance, 

PsourceC'^GLli ^os) X pions (^GLi i J-'ol) ■ Contour levels are separated by factors of 10, the dashed contours mark a density of = 10^* (p 
in units of solar masses per cubic parsec). 

in the density contours arise from the maxima in the source density and lens density, respectively. It can easily seen 
that lowering Dqs would slightly decrease the needed lens-source separations, but at the same time the source density 
decreases rapidly. 

For bulge sources at 770 kpc in Fig. [7] a,b the minimum source-lens separation is Dos — Do\ = 74 kpc for a low 
mass lens with Mq = O.IM©, 16kpc for a O.SM© lens, and 8.1 kpc for \Mq. For masses lower than O.SM© the large 
lens-source distance required to create the observed lensing event, implies a very small event probability, because for 
large lens-source LOS distances either the source or the lens density is small. So, from just using the brightness of the 
event, lensing of bulge stars < O.SM© is very unlikely. For a bulge source more reasonable lenses for self lensing are 
high mass bulge remnants (Z?os — Do\ = 0.16 kpc for 50Mq) or very rare, young, high mass disk stars (0.82 kpc for 
IQM0, 1.6 kpc for 5Mq). 

Despite the lower disk density at the position of GLl disk sources lensed by bulge remnants or high mass disk 
stars seem more plausible (see the two right panels of Figs. [T]: and[7ti). We used a reasonable disk extinction of 
Ar = 0.68 mag with the source color and brightness derived from the light curve fit. These values are only consistent 
with a population older than > 0.5 Gyr (see blue marker in Fig. [2k)- Therefore for these disk sources the maximum 
Fo/R^-T&tio (at {R — I) = 0.59 mag) is obtained by stars with Mr — —1.83 mag and i?* = 35. li?©. 

For a disk source at 770 kpc a lens with IMq has to be separated by 3.7 kpc; for a lens with Mq = ^Mq this 
separation becomes 0.75 kpc. The distances are smaller, since the disk sources are intrinsically brighter, and thus a 
lower magnification is needed for the event, and therefore, lenses can be closer to the source. 

Of course, the event rate contribution is not just proportional to the mass density product of all lenses and all 
sources (the product of that is displayed as contours in Fig.l?]) but depends on the number density of those sources and 
lenses only that can produce the observed event. There, the (M/L), and other stellar population properties (fluxes of 
post-main sequence stars, fraction of stars in a certain color interval) enter. 

However our qualitative result, that WeCAPP GLl lenses below O.5A/0 for bulge sources and below O.IMq for disk 
sources are very unlikely, is consistent to the quantitative result in Fig. [8l 

5. QUANTITATIVE WECAPP-GLl ANALYSIS INCLUDING FINITE STELLAR SIZES 

In ^J3]we treated stars as point sources, like iRiffeser et all (l2003f). We n ow abandon the point source appr oximation 
and use for each star the radius based on the isochrones from iGirardi et al. (2002). We evaluate Eq. lCTll using the 
stellar source size distributions, and show the results in Fig. [S] 

Figure [8] shows the halo lensing and self lensing event rates (per flux excess and FWHM-time of the event and 
per year and square arc minute) at the location of WeCAPP-GLl using its brightness, time-scale and color and the 
extended stellar sizes. Comparison to Fig. |4] shows, that self lensing becomes dramatically suppressed; the vertical 
scale changed by one order of magnitude relative to Fig. H) The halo-bulge lensing rate decreases by a factor of 6 
relative to the point-source formalism. The self lensing rate however decreases much more dramatically using the finite 
source description. Relative to the most likely M31 halo-bulge lensing, self lensing (bulge-bulge, disk-bulge, disk-disk 
and bulge-disk) is about 2 orders of magnitude less likely, e.g. lensing of a bulge star by a O.8M0 MACHO in M31 
(black points) is 119 times more likely as bulge-bulge (red curve) or bulge-disk [green curve) self lensing. The total halo- 
lensing contribution (summing over M31 and MW halo-bulge and halo-disk lensing) is 63 times more likely than all self 
lensing contributions assuming a 100% MACHO halo. In addition, for each lensing configuration, the probability for 
small lens masses is reduced. This can be explained as follows: The large magnification needed for bright events can be 
obtained in the point source approximation either by large Einstein radii (efficient source-lens geometry or large lens 
masses) or events with small impact parameter. For point sources, the magnification can formally become infinitely 
large if the source passes the LOS to the lens. For extended sources the magnification saturates once the source comes 
closer to the lens LOS than its projected source radius. The source radius projected onto the lens plane is the larger, 
the closer the lens-source pair, and therefore self lensing events are most strongly suppressed in magnification and thus 
in flux excess at the light curve maximum. The only way to obtain bright events with finite source sizes is to have 
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Fig. 8. — Event rate r{Af) (per area, per event brightness A_F, per event time scale tpwHlvii per color (R — I)) for WeCAPP-GLl using the 
evolved MF including remnants. The lens mass prolaability was determined allowing for finite stellar sizes, using the measured ip^fji^i) 
^mcas g^j-jjj ^moas q£ ^[^g event [{R — I)o = 0.59 for disk sources, {R — I)o = 0.70 for bulge sources, with a Gaussian error of 0.1], assuming 
a mean M31 extinction of 0.19 and 0.51 for bulge and disk stars, and a MW extinction of 0.17 (in the R-band). The red, green, blue and 
purple curves show the lens mass distribution function for bulge-bulge, bulge-disk (bulge lens and disk source), disk-bulge and disk— disk 
lensing. The information about halo lensing (M31-halo-bulge and M31-halo— disk data points are in black and grey, and MW-halo-bulge 
and MW-halo— disk in brown and orange) is displayed as points and not by a curve: we assume a mono-mass-spectrum for the halos, and 
each point represents a halo which is made of a given MACHO-mass to 100%. 

spatially well separated source and lens stars (increasing the Einstein radii and decreasing the projected stellar sizes) 
or to have large lens masses (increase of Einstein radii). Therefore O.IMq (for M31-halo-bulge lensing) and COIM© 
(M31-halo-disk lensing) MACHOs have too small Einstein radii to achieve the required large magnification once finite 
source sizes are taken into account. The suppression sets in for smaller lens masses in the M31-halo-disk lensing case, 
because disk stars are brighter than bulge stars. 

To illustrate the differences between point source and extended source sizes approximation in more detail, we evaluate 
Eq. (jCSP and separate the events into those which do not show any finite source signatures in their light curves (no 
fss) and those which show finite source signatures (with fss) and display the corresponding results in Figs. [H and [TOl ^ 
To understand the differences to Fig. [5] it is useful to recall the relations for event time scale tpwHM and magnification 
of events at light curve maximum Aq in the high magnification approximation [for carrying out the integrals we used 
the exact relations described in lRiffeser et al\ (j2006[ )]: 



^fwhmI^) = Whm(&) oc tEUQ = tE^ = vt , b>b* 

tpWli! (b) = ipWHM ib) ~ tFWHM (b) ^ Vi2b*/by - 1 > tpWHM (b) , b<b* 

f-(&) = Ao(6) « i?E/fo , b>b* 
A^'^'"^''^{b)=Ao{b*)KRE/b* , b<b* 



(5) 



(6) 



with b* ~ i?*I?oi/(2-Dos)- In the point source approximation, source lens configurations with small transverse velocity 
and small impact parameter can contribute to the event rate at given tpwHM; in the finite source treatment, transverse 
velocities have to be larger than a limit cx 6*/iFWHM to make events without finite source signatures (no fss); in 
the other case, Vt < , light curves will show finite source signatures (with fss) (i.e. show a saturated magnification, 
and an increased event time scale). Therefore, Fig. [TO] can be constructed out of Fig. O by moving the event rate 
contributions due to events with vt < to events with finite source signatures, after accounting for their increased 
time scale, and their decreased brightness. By comparing the panels in the last row in figures Figs. [5] and [9] one can see 
that finite source effects are almost unimportant for Milky Way MACHO-events with lens masses between 0.2Mq and 
IOOMq: for both MACHO masses, the event rate and the characteristics of the events with time scales larger than 1 
d (which are observable with present experiments) are nearly the same in the point source approximation and in the 
extended source size treatment. The remaining panels in Fig. [9] show, that finite source effects are less important for 
heavy M31-MACHOs, and are relevant for low mass M31-MACHOs, like 0.2Mq (second row sub panels), and that 
they are extremely important for the correct interpretation of the self lensing contribution (first row sub panels): an 
event like WeCAPP-GLl can hardly be caused by MACHOs with masses much smaller than O.2Af0, otherwise the 
flux-excess cutoff would be lower than the event brightness, and WeCAPP-GLl is unlikely to be caused by self lensing. 

To convert to the same units as in Fig. [8] use Eq. IC8I E.g. in Figs. l9l and [TOl the WeCAPP-GLl measurements cross the contour 
level of haloM3i(A^=0.2A^Q)-bulge lensing at -2.30 and -1.86, respectively. This converts to r/[yr~^ arcmin"^ mag"'^] = (10"^ '^'' + 
10'i *'5)/{0.41n(10)2 ■ 1.83day • 10.07 X 10"^ Jy) X 10"-* = 4.8 X 10"^ (lO""* Jy)-^ day-^. 
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Fig. 9. — Flux excess and time scale distribution of events with WeCAPP-GLl color {{R — I)o = 0.5 9 for disk sources, [R — I)o = 0.70 
for bulge sources, without using its error) and its location. We make use of finite stellar sizes in Eg. ICSI but display only events which 
do not show any finite source signatures (no fss) in their light curves (left part of brackets in rhs part of Eq. The two further 
obscrvablcs of WeCAPP-GLl (flux excess and event ti me scale obtained fr om a point source lensing light curve fit) are marked in green. 
The estimate for the WeCAPP 6a detection limit (see IRiffeser et oil I I2006I ). Table 1) at the position of GLl is shown as white line. The 
contours differ by 10""'^, and the dashed contour is that of 10~^ events per year, per square arc minute at location of WeCAPP-GLl, per 
color (in magnitudes), per logarithmic timescale (time scale in days) and per flux-excess (converted to magnitudes) (see footnotes [3] and [Jjl . 

Among all self lensing scenarios (Fig. [9l first row) the brightness of WeCAPP-GLl is most easily to achieve by bulge 
lenses (Fig. [SI first two panels of the first row), because in that case source and lens stars are relatively separated (more 
than for disk lenses). 
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self lensing 
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0.0101 
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< 19.0 


0.0204 


49.0 




> 0.9 


< 21.3 


6.47 


0.155 


halo lensing with 20% IMq 


> 10 


< 18.7 


0.0541 


18.5 




> 8 


< 19.0 


0.0969 


10.3 




> 0.9 


< 21.3 


2.21 


0.45 



TABLE 3 

Expected rates for short time scale (FWHM-timescales between 1 and 3 days) lensing events in the entire 17.2' x 17.2' 

WeCAPP, assuming 100% detection efficiency. 

In Table 3 we have evaluated the number of events we expect to find with time scales between 1 and 3 days in 
the entire 17.2' x 17.2' field monitored by WeCAPP (assuming 100% detection efficiency). Short self lensing events 
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Fig. 10. — Flux excess and time scale distribution of events with WeCAPP-GLl color [{R — I)o = .59 f or disk sources, (i? — /)o = 0.70 
for bulge sources, without using its error] and its location. We account for finite stellar sizes in Eg. ICSI but add only events (right part 
of brackets in rhs part of Eq. which do show finite source signatures (with fss). The two further observables of WeCAPP-GLl 
(flux excess and event time scale obtained from a point source lensing light curve fit) are marked in green. The estimate for the WeCAPP 
6(7 detection limit (see IRiffeser et oil I I2006I) . Table 1) at the position of GLl is shown as white line. The contours differ by 10""'^, and 
the dashed contour is that of 10~^ events per year, per square arc minute at location of WeCAPP-GLl, per color (in magnitudes), per 
logarithmic timescale (time scale in days) and per flux-excess (converted to magnitudes) (see footnotes \3\ and [2J ■ 

with a peak magnitude^ larger than 19.0 mag (as lower limit for WeCAPP-GLl, see Q iii the should take place only 
every 49 years. For reasonable halo mass fraction of 20% consisting of IM0-MACHOS every 10 yrs (short) events 
above a 19.0 mag threshold could be observed. These values decrease very rapidly lowering the threshold to the most 
likely peak brightness for GLl of 18.7 mag: The mean time between events for selflcnsing is 99 yr whereas for halo 
lensing 18 yrs. While for the short, fainter events (< 21.1 mag) self-lensing would dominate by a factor of 3. Note 
that these values represent the theoretical event rates in the WeCAPP field without detection bias. In other words 
only if a survey is complete to a certain threshold it should be possible to measure these numbers of events. For real 
microlensing surveys the detection efficiency for bright events will be much higher than for faint events, and therefore 
the ratio of bright vs. short events would be different. If we assume to have the same detection efficiency for halo and 
selflensing (which in principle can differ according to there intrinsic distribution) the difference between halo-lensing 
and self lensing favors the assumption that WeCAPP-GLl {green line) was caused by a halo lens. 

To summarize, an event as bright as WeCAPP-GLl is extremely unlikely to be a self lensing event. The most likely 
self lensing scenario is bulge-disk self lensing. MACHO-lenses above 0.2 solar masses are expected to cause the bright 
events much easier. 



6. INTERPRETING WECAPP-GLl: IMPACT OF EXTINCTION 

IRiffeser et al\ (|2006l ). which assumes an on average M31-extinction for all disk stars in the WeCAPP field of 0.51 mag 
and for all bulge stars of 0.19 mag in the R-band, independent of the angular position of the event, and independent of 
the LOS distance to the source. The MW-extinction is set to 0.17mag (R-band). The result for this M31 standard-dust 
model has been described in Fig. [5] already. 

To see to which extent dust can change the interpretation of lensing events, we set the M31-extinction (not the 
MW-extinction) along the LOS to WeCAPP-GLl to zero. We do not assume that the mean M31-extinction of bulge 
or disk stars goes to zero as well^, but that by accident just this one sight line to the source is not extincted by 
M31-dust at all; the extinction then equals that of the MW of about 0.17 mag. The results are shown in the left-most 
column of panels in Fig.[TlJ With no M31-dust along the LOS to WeCAPP-GLl the event would become intrinsically 
fainter (relative to the extinction case and giving the observed fluxes) by 0.19 mag and 0.5 magi magnitudes in the 
R-band for bulge and disk sources. This reduces the necessary magnification and makes self lensing more likely. At the 

^ We define the function magj5(x) to replace the transformation from fiuxes to magnitudes magjj^(x) = — 2.51og]^Q(x/i<Vcga,R) with 
i^Voga.R = 3060 Jy. The inverse is defined as fluxR(a;) = i^Voga.RlO"" *^ . 

This would imply a decrease of the disk and bulge number density and thus of the event rates for disk source stars by 37% and for 
bulge source stars by 16% and thus would decrease the lensing rate of disk sources relative to bulge sources 
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same time the lensed sources are expected to be intrinsically redder, which changes the stellar types of stars and thus 
its absolute brightness {lower panels). One can derive the increase of the event rates for an event like WeCAPP-GLl 
also by shifting the contours in Figs. [51 and [TUl by 0.19 mag and 0.51 mag along the Ai^-direction for lensing of bulge 
sources and disk sources respectively. The ratio for the event rates of all self-lensing to all halo-lensing configurations 
is still of the order of 1:100 for full MACHO halos and MACHO masses in the range of 0.1 to 1 solar masses. So even 
in the implausible case of no M31-dust along the LOS to the source of WeCAPP-GLl, a 1% MACHO fraction of the 
M31 and MW halos already provides the same lensing rates as self lensing does. 

As alternative values for the total (MW-I-M31) line of sight extinctions we also use values of 0.7 mag and 1.05 mag 
("strong extinction") and show the results in the third and fourth columns of Fig. [TlJ In this case lensing rates are 
suppressed and the most likely lens masses are shifted to higher masses. This can be understood since now only higher 
mass lenses are able to provide the high flux excess found for WeCAPP-GLl. 
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Fig. 11. — Event rate r(A/) (per area, per event brightness AF, per event time scale ipwHMi per color (R — I)) for WeCAPP-GLl. 
upper panels: without changing the color by the assumed extinction [{R — I)o = 0.59 for disk sources, (R — I)o = 0.70 for bulge sources, 
with a Gaussian error of 0.03], lower panels: taking into account the color shift due to the assumed extinction. The red, green, blue and 
purple curves show the event rate for bulge-bulge, bulge-disk (bulge lens and disk source), disk-bulge and disk-disk lensing. The information 
about halo lensing (M31-halo-bulge and M31-halo-disk data points are in black and grey, and MW-halo-bulge and MW-halo-disk in brown 
and orange) is displayed as points and not by a curve: we assume a mono-mass-spectrum for the halos, and each point represents a halo 
which is made of a given MACHO-mass to 100%. The results have been obtained with an average dust extinction of total bulge and disk 
population of ytjj = 0.36 mag and J^n = 0.68 mag. These values are expressed in the two middle panels as 0.35 and 0.7. The results for 
the minimum line of sight extinction of 0.17 (set by the MW) and for a high extinction of 1.05 are shown to the left and right of the same 
row. The 4 panels in the first row show results analogous to the 2""^ row, with the difference, that the de-reddening of the source star has 
been neglected: this simplification would have a minor impact on the interpretation for the stellar populations assumed here. 



In the second row of Fig. [TT] we show how the extinction influences the color of the event and therefore the type of 
source stars which were most likely lensed. Whereas in the upper panels the event rates are shown, if one accounts 
for the dimming but not for the reddening by dust, in the lower panels we take into account the reddening with the 
assumed extinction. The effect is lowest for low extinction; for high extinctions the predicted masses are strongly 
changed. The assumed shifts in color {R — I) for event extinctions of 0.17, 0.35, 0.70, 1.05 mag are -0.06, -0.12, -0.25, 
-0.37, respectively, which means that the intrinsic colors of GLl are 0.77, 0.71, 0.58, 0.46. Comparing these colors 
with the CMD in Fig. [2l shows that the intrinsic brightness can drastically change for slight color differences: for 
the 12.6 Gyr component from -1.0 to -1-3.7 mag, for the 2.5 Gyr disk component from -1.8 to -1-2.3, whereas the 20 
Myr disk component slightly brightens from -6.2 to -6.9 mag. This can explain the extreme differences in the lens 
mass distributions and underlines the importance of accurate color measurements and precise extinction estimates. 
Although the color plays an important role in the correct mass interpretation, the ratio between halo- and self-lensing 
is almost independent of the assumed extinction. For all extinction assumptions WeCAPP-GLl is more likely caused 
by a halo lens than by a stellar lens. Only for the de-reddened, high extinction case {last panel, 2"'^ row) the halo to 
self-lensing ratio approaches a factor of 10. 
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7. THE IMPACT OF A METAL POORER DISK 

The maximum flux excess AF of an event depends on the flux and radius of the source star. Therefore, the question 
arises, how sensitive the interpretation, that WeCAPP-GLl can hardly be achieved by self lensing, depends on the 
assumed metalicity of the stellar population. The luminosity and sizes of post-main-sequence stars change slightly, if 
the metalicity of the composite stellar population is altered. 

We therefore have modeled a metal poorer disk (|Williamsl |2002[ ) by changing the metalicities to Z ^ 0.008 (see 
Girardi's isoc_z008.dat). The change of metalicity has several competing effects: 

i) a metal poorer population is brighter 

ii) for a metal poorer population, smaller MACHOs masses are allowed, since brighter stars need smaller magnifi- 
cation, which can be produced also by smaller masses without finite-source-size saturation of the magnification. 

iii) the lens properties are nearly not affected by the luminosity evolution, since reducing the metalicity mostly 
changes the properties for the lensed stars (sources), and mass evolution takes place only for very few stars. 

iv) a stellar population with lower metalicity (changed from solar by a factor of 2) contains PMS star brighter by 
a factor of roughly 2 (« 0.75 mag). Their radii (i?* c>c \/Fo) are larger by a factor of roug hly a/2 = 1.4. bmce 
(see Eq. (j31)) AFinax oc Fo/i?* oc ^/Fq oc R^. for events with the same color, one expects the flux excess of the 
brightest events to increase by a factor of 1.4, or 0.38 mag if the stars become metal poorer. This would make 
the contours in Figs. [H] and [TO] to shift by 0.38 mag in vertical direction and push the flux-excess limit due to 
finite source effects by the same amount. 

v) However, if stars are brighter, fewer of them are needed to account for the total observed light, which reduces the 
overall lensing rate, e.g. the (M/L) of the disk population drops from 1.2 to 1.0 in the R-band if the metalicity 
of the population is changed from solar to 0.008. 

All these points explain the relative small differences in the event rates by changing the metalicity of stellar popula- 
tion. However the lower part of mass probability distribution for a measured event can slightly shift to lower masses, 
as the lower amplifications can be produced by lower masses. 

8. OUTLOOK: STATISTICAL INTERPRETATION - HALO VS. SELF LENSING FROM BRIGHTNESS DISTRIBUTIONS OF 

EVENTS 

We now show that the brightness and time scale distribution of events can be used to discriminate halo and self- 
lensing. 

Figure [12] {left panels) compares the event rates as a function of event brightness for self lensing and halo lensing 
events within rings around the M31 center. These rings have a thickness of 1 arcmin and outer radii between 1 arcmin 
and 12 arcmin (blue/red for self lensing, green/black for halo lensing). The event rates were calculated ignoring 
detection efficiency factors, assuming a 12.6 Gyr bulge and a composite disk population, and the mean extinctions for 
the bulge and the disk sources. We consider only events longer than 1 day (ipwHivi)- For self lensing, the event rate 
drops steeper with event brightness within central rings than within outer rings. This is caused because in the inner 
rings self lensing of bright events is suppressed by the closer LOS distances of source-lens pairs. For the outer rings the 
mean distances between disk-bulge pairs are larger and the lensing rate is increased. For halo lensing there is no such 
effect. The brightness cut-off of the event-rate-event-brightness relation is determined mostly by stellar population 
properties (self lensing) and by the MACHO mass (halo lensing), as can be seen by comparing the left panels in 
Fig. [T^ which were obtained for O.IA/q, 1.0A/q, and IOA/q MACHOs from top to bottom. This demonstrates that 
the brightness distribution of bright events could itself discriminate halo and self lensing if a sufficient large number 
of events is available. The two lines represent the WeCAPP-GLl flux excess (fit and highest data point). 

The distribution of event time scales (only events brighter than 10~^ Jy in R are considered) is fairly similar for self 
lensing and MACHO events (in slope) if the MACHO mass is O.lAf© (this can also been seen from Fig. [9] which shows 
for the location of GLl that the event characteristics [flux excess and time scale] are most similar for self lensing and 
halo lensing if the halo population has a mass of O.2M0). However, if MACHO masses are larger, events will have 
longer timescales than predicted from self lensing (see Fig. [T2l right panels). 

The right panels in the 2""^ and 3"''' row of Fig. [TO] show the time-scale distribution for IMq and IOMq MACHOs, 
respectively. Because of the clear difference between selflensing and halolensing the time scale-distribution is a good 
discriminator if the MACHO mass is larger than l.OAf0. 

The sharp cut-off in the brightness distribution of halo-lensing events is caused by the fact that the MACHO mass 
function used is a deZta-function. Since the different lensing masses produce a brightness distribution which has quite 
a different gradient with respect to selflensing, it is conceivable that there exists a specific halo-mass function where 
the brightness distributions look similar to self-lcnsing. Therefore, one probably can find a halo mass function that 
can mimic self lensing as long as only the brightness distribution of events is considered. If one however combines 
the brightness distribution with the distribution of the time scales and locations of the events, the halo-self-lensing 
degeneracy can be broken. 



10+1 
10+0 
10-1 
T 10-2 
1 10-3 
10-'' 
10-E 
10-6 
10-^ t 



10- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

rHalolensing: 0.1 l^sun 
rSelf lensing 



WeCAPP-GLl: evidence for a MACHO component in M31? 

10+"! 



15 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



eElLii 




R=o-r ^ 

Ffc10-11" 

ii III I fl=1 1-121 
I II ill II li Ill I I I r 



10+1 
10+« 
10-1 
10-2 
10-3 
10-1 



17.5 18.0 18.5 19.0 19.5 20.0 20.5 21.0 
-2.5 log &F 1 3060 Jy 

l|ii II II II i|i II II II ii| 



laioifeiwihy; 1 .0 
^Self lensing 



II III |i II II III i| II II III II |i II II III i| II II III ii 





I Ii I I Inn 



FUO-r ^ 
ti -2' 

I 

^4-5" 
tE-e' -! 

tiO-11 _ 
Ffci1-12'" 

I r 



17.5 



18.0 18.5 19.0 19.5 20.0 20.5 21.0 
-2.5 log AF / 3060 Jy 





10 




10 


> 


10 


a) 


10 




cc 






10 




10 




10 




10 



10+1 

10+« 
1 

2 
3 
4 
£ 



II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
sun 



^^Haiolensing: 10 
^Self lensing 




Ffco-r ^ 



^4-5" = 
tS-S' -s 

1=7-3' -. 

■fl RI I I II li Ill I I I r 



17.5 18.0 18.5 19.0 19.5 20.0 20.5 21.0 
-2.5 log aF / 3060 Jy 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

H=0-1' : 




10 15 20 25 

t^Pihm[c^] 



30 



10^1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

R=0-1' ; 




10 15 20 25 

tfwhm [ d ] 



I I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

R=0-1' : 




10 15 20 25 

ttmhra [ d ] 



30 



Fig. 12. — Left panels: Brightness distribution of events (with tpwHM > Irf) within rings of 1 arcmin around the M31 center. Outer radii 
are chosen between 1 arcmin {red for self lensing, black for M31 halo lensing) and 12 arcmin {blue for self lensing, green for halo lensing). 
The two vertical lines {grey) represent the GLl flux excess (fit and highest data point). Right panels: Time scale distribution of events 
(brighter than AF > 10"^ Jy corresponding to an excess magnitude of 21.2) within rings of 1 arcmin around the M31 center. Outer radii 
are chosen between 1 arcmin {red for self lensing, black for M31 halo lensing) and 12 arcmin {blue for self lensing, green for halo lensing) 
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9. DISCUSSION AND CONCLUSION 



We showed that accounting for extended sources can dramatically change the lensing and self lensing rates for events 
as bright as WeCAPP-GLl. The reason is that magnification saturates, and that the maximal brightness of an event 
for a given source size and source flux depends only on the source and lens distances and lens mass (not on the impact 
parameter): Very bright events thus require a minimum source lens distance, which makes halo lensing more likely 
relative to self lensing. For WeCAPP-GLl the inclusion of the finite stellar sizes makes the ratio of self lensing to halo 
lensing by about a factor of 8 less likely compared to the point source approximation. This result was obtained based 
on simple stellar population models of the bulge and disk and on a simple description of the dust extinction and is 
almost independent of MACHO mass as long as that is in a range of 0.1 — 3.5Mq. 

Likewise if one uses position, FWHM time scale, flux excess and color of WeCAPP-GLl, self-lensing is even 13 times 
less likely than lensing by a MACHO, if 1 solar mass MACHOs contribute only 20% to the total halo mass in this 
mass range. 

Assuming a 100% detection efficiency we expect a self lensing event with with FWHM time-scales between 1 and 3 
days and a peak flux excess brighter than 19.0 mag in the entire WeCAPP survey field only once every 49 yrs (99 yrs 
for magf{(AF) < 18.7 mag). On the other hand if 20% of the halos of M31 and the MW are made of 1.0 solar mass 
MACHOs then an event like WeCAPP-GLl would occur every 10 years (every 18 yrs for a 18.7 mag threshold). 

This implies that a small fraction of baryonic MACHOs like (brown dwarfs, neutron stars, stellar black holes, cool 
white dwarfs) in the M31 halo is sufficient to make WeCAPP-GLl more likely to be a halo-lens event. 

We have also shown that different values for the extinction of the WeCAPP-GLl event do only slightly change the 
halo-lensing - self lensing ratio, but more strongly change the most likely lensing masses. This is because the major 
impact of the extinction is that of a bluer intrinsic source color and therefore implies a change of stellar source size 
and brightness. 

We emphasize that any interpretation for high difference flux event s and their predicted rates has to account for 
finite-source effects. Further published examples are lAnsari et al\ (|1999f ): lAuriere et al\ ()2001| ): [Paulin-Henriksson et al\ 
([2003). Therefore the brightness distribution of events in general is a good discriminator between self and halo lensing. 
The time-scale distribution is a good discriminator if the MACHO mass is larger than I.OMq. 

Our model reasonably describes the {M/L) for the stellar population in M31 and we are confident that the mass 
density and velocity distribution of the stars is fairly well modeled. With these ingredients one cannot obtain an 
event rate high enough to make a detection of WeCAPP-GLl likely within our survey interval. Therefore, our analysis 
supports the existence of a (small baryonic or non-baryonic) MACHO component in the halos of M31 or MW. 

In this paper we have ignored detection efficiencies. We were allowed to do so, because, for a single event, the 
probability ratio of halo lensing to self-lensing is independent of detection efficiencies. For WeCAPP-GLl, halo lensing 
(assuming a 20% MACHO halo) is 5 times more likely than self-lensing. We admit that this factor of 5 is not large 
enough to conclude that self-lensing is definitely excluded. However, we point out, that 'WeCAPP-GLl '-like self-lensing 
events are very rare. Put in another way, the fact that this event was observed can be more easily understood, if there 
is a, e.g., 20% compact dark matter fraction in the halos of M31 and the MW. The analysis of more bright events, and 
the comparison of the event rates for bright and faint events will give more insight, as well as the use of their time scale 
distribution. This however requires the knowledge of the detection efficiency (as a function of brightness, timescale, 
location and color) of the events, and will be subject of a further paper. Our results provide a strong motivation to 
search for the brightest short time scale events towards M31 in wide field surveys (like Pan-STARRS 1). 



Assuming an M31-distance of 770 kpc we converted our light model (described in lRiffeser et al\ ()2006D ) into mass 
using the theoretical bulge and disk R-band (M/L)]^ of 3.3 and 1.2 derived using their theoretical stellar populations. 
The colors extracted from these stellar populations are for the bulge {B — R)q — 1.63 [{B — i?)meas = 1-91] , {V — R)o — 
0.60 [(y-i?)„,oas = 0.72] and for the disk {B - R)o - 0.81 [(B-i?)„,eas = 1-33], {V - R)a = 0.37 [(^-i?)n.cas = 0.60]. 
The analytical mass density models are 



APPENDIX 



M31-MODEL 




a < 3.1 pc 

3.1 pc < a < 20 pc 

a > 20pc 



(Al) 



with a EE 0.57^2 + y/0.5T^z^ + .t^ + y'^ + l.llz^ in kpc, 



Pdisk{x,y, z) = O.27Af0/pc^ x exp 



( 




(A2) 



with her = 6.4 kpc and hz = 0.3 kpc. 



Phedoix, y,z) = O.O65M0/pc^ x 



1 



r < 100 kpc 



(A3) 



1 + {r/r,f 
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^'^.(-t) = z^-t exp H^], (A4) 



with r = {x"^ + y"^ + z^)^!'^ and = 4kpc. 

Integrating over these density profiles gives a total mass for the bulge of 4.44 x lO^^M©, for the disk of 4.18 x lO^^M© 
and for the halo of 122.7 x VS^^Mq. 

We used a bulge and disk incHnation of i = 77°, a position angle of the disk major axis of P. A. — 38° and P. A. — 50° 
for the bulge. 

For the transversal velocity distribution we assumed 

1 / ^vl\ (vQVt 

— vt exp — — MO ~ 

with (Tis — [(Ti^ + (cTgDoiZ-Dos)^]"'^- The velocity dispersions are cTbuigc ~ 100km s^^, Cdisk = 30km s~^, 
Chaio — 166 km s^ , cTMW-haio = 156 km s^ . The additional rotations are taken into account as 
uo(-Dos, Z^oi, t^rot,i, I'rot.s, Wo-M3i), whcrc we used Wiot, bulge = 30 km s~^ and Wrot,disk = 235 km s"^ and an observer's 
motion of Vo-msi = 129km s"^(Dos - -Doi)/I?os- 

EVENT RATE 

Eq. [5] contains the following analytical functions [we have dropped the variables on the right-hand-side (rhs) mostly; 
all functions on this side can be expressed as a function of the variables on the left-hand-side (Ihs)] 

T(AF/^o + l) = T(Ao) = 2^u{^)' - .(^o)^ = V8 "^tZ^^^^^^^^^': 

V|/(Ai'/^0 + Ij - JL (Ao)-4V2 (Ao^-l)V4(A„+3)i/2 

D:,{R.,D,.,M,^F,Fo,) = D,,(l + ^^^fl±^y' with 
u*o{R,,DohDos,M) = 




n*iAF,Dos,Fo,M,R,) = ^^^^p^ =2CDos'{Fo + AF){CDos + AFi2Fo + AF))-^ 



T*{u,,R.,D,,,D,^,M) ^ 2j-j^^g^^~2~u,^ with A5 = ^l+(2|^ 

i;t(i^oi,Cos,M,iFWHM,Ai^,Fo) = ^g;';;-^ T 

<(uo,£'oI,£'os,M,tFWHM,i?*) = J^T*{UQ,..) 

and the color-magnitude relations of stars, Pcmd(A^,C), the mass function of stars and potential MACHOs, S,{M), the 
spatial density of sources n{x, y, Dos) on the line of sight, the mass density of lenses, p{x, y, Doi), and the transversal 
projected lens-source velocity distribution Pvtivt). The mass density of lenses depends on the dynamical model and 
the number density of sources is constrained by the observed light of the stellar population. If the extinction value of 
the stellar population changes, or equivalently, the {M /L) of the stellar population is changed, the event rate scales 
linearly (and not quadratically). 

The first term inside the brackets of the rhs of Eq. [5] collects contributions from events that do not show finite source 
signatures, i.e. ev ents with impact par ameters larger than the projected source radius, uq> Uq~ R^,Do\ {2R^Dos) (see 
Eqs. 65 and 66 in lRiffeser et al\ ()2006D V The second term on the rhs is due to events showing finite sou rce signatures, 
i.e. ev ents with impact parameters closer than the projected source radius, uo < Mq; used Eq. 67 in lRiffeser et al\ 
(HH). 

All functions with depend on the radius i?*(A^,C). 

LENS MASS PROBABILITY DISTRIBUTION 

The event rate differential in Eq. ^ depends on 7 arguments: The first 5 arguments (a;, y, ^fwhMi C) are 
observables of lensing events; the brightness of the source stars, A1, is not directly observable in light curves (the color- 
luminosity distributions of stars can instead be statistically described using observations of resolved stellar populations 
in M31 or using stellar population models); the sixth argument, the lens mass, M, is the quantity to be statistically 
constrained (lens mass functions and amplitudes) from observing lensing events. One can integrate Eq. ([2]) to obtain 
lower order differentials: 

d'^T{x,v,tF^umAF,C,M) ^ f d^r(x,y, Whm,AF,C,M,A^) , . . (CI) 

dx dy dtpv/HM dAF dC dM J dx dy dt-ev/um dAF dC dM dM 

d°r(a:,y,tFWHM,AF,C) _ r r rf^r(z,t/,tFWHM:AF,C,Af ,A^) ^iK AJn^ (C2) 

dxdydtj^vjiiMdAFdC ~ J J dx dy dt^viUM dAF dC dM dM ""'^ 
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dx dy dtpwHM dAF ~ J J J dx dy dtpwHM rfAF dC dM dM "-^ ' 

In the following we assume, that one can measure the location of an event without any error. Some of the remaining 
observables (^fwhivi, Ai^,C) might have fairly large errors. Let p{o, o™'''*^) be the probability for measuring 0"^°^^ with 
o being the true value. We then can estimate the error-weighted values of the differentials in Eqs. ICl[IC2l and lC3l at 
the location of the observables: 

• without color measurement 



' d^r{x,y,tFWiiMAP) \ (4- Ap\ ^;4p(^mca, ^^moas ^tpyf^MAF) ^^^^ 

dxdydtpvjHM dAF /^,„<„3 „moas l''FWHM, ) ■ dxdydtp^uM dAF 



I rf^r(3:,y,tFWHM,A_F) \ ._ 

\ dKdydtpwHMdAF / Apmcaa'" ' 

' X ,y ,i,p^jj^,tij' 

/ / ''"^C;'CwHM^AF^'^^ p(ipwHM, ^^^^Im) P(A^^, A—) dtpwHM d^F 



rf^r(a;,y,tFWHM,A_F,M) \ ( A J\ i np.\ 

dx dy dtpwuM dAF dM , „ ) — I'-^OJ 



///- 



r(x'"°°'°,i/'"°'^°,tFWHM.AF,M) j-mcas \ n( \ T? A™oasN Jj- rl A P 

dxdydt^y^sM dAFdM Pl^FWHM , IpWEM ) P(^^ ^ ^F ) aiPWHM a/\-t 



• with color measurement C™'^'^'^ 



d r(tFWHM,AF) (C^7\ 

d^r{x,y,tpwHM,AF,C) \ /, \'P\ •— r d°r(a:°'°''°,i/'""'°,tFWHM,A-F,C) /p naicas^ jn 

dxdydtFWHMdA_F<ic/^,„,,, ,^„„as^C"-™ ™ dxdydt^^amdAFdC f^^^^ > 

Eg. lCTI can be transformed (To transform a distribution df{x)/dx to dg{u)/du with u(x) and its inverse we 
used dg{u)/du = (i/[a;(w)]/(ix |(ia;(w)/(iu|, see also footnote [3|) into 

^|lj«=„.4 1„aO,MO»flu.„,«^£&5«^ ,C8, 

dOm UlFWHM aAi* 

defining = logio^pwHM and J„ = magj,(AF). 



d^r(3;,y,tFWHM,AF,C) \ 

da; dy dtFWHM dAF dC / „e„ „moaB rmoas f mcas Amcas' \ ' 

/ / / '^'^drdydCwHMlAFdy-'^ ^'(iFWHM, t^^^l^yd p{AF , A-<=-) p(C, C"--) dtpwHM dAF dC 



d r(a:,y,tFWHM.AF,C,Af) \ r ._ /'Plfl'l 

da; dy dtFWHM dAF dC dM I, , . „eas ^„oas ^ ' ^'^ ) 

in '^^^TxTyfZ.ll^lFdfdM''''' P(iFWHM, t^^\^) p{AF ^ A^^) p(C, C^^) dipwHM dAF dC 

Equation (jClOp gives the contribution to the event rate (with event characteristics as observed) as a function of lens 
mass. The probability for a lens with mass M causing an observed event can therefore be written as 

\ dx dy aipWHM dlXb dC dM / „eas ,,meaB ,moa= A^^^., f,„,oaB 

We will evaluate the lens mass probability function for each self lensing (bulge-bulge, bulge-disk, disk-bulge and 
disk-disk) a nd h alo (M3 1-halo-bulgc, M31-halo-disk, MW-halo-bulge, MW-halo-disk) lensing scenario. 

Equations IC7I and IC4I describe the distributions of events (at the location of the observed event) in the flux-excess- 
lensing time-scale plane, for the case where the observed color of the event is used (Eq. IC7|) and for all events, 
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irrespective of their color (Eg. IC4[) . Studying the event distribution in this plane is very useful, since one can immedi- 
ately see, if an event is unlikely (given a theoretical model), and how data quality (imposing limits on detectable flux 
excess and time scale) restricts the measurable event rate. 

The relative values of Egs. IC9I and IC5I give the relative probabilities for the different lensing scenarios, if the color 
of the event is used or not used, respectively. 
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by the Sonderforschungsbereich SFB 375 of the Deutsche Forschungsgemeinschaft (DFG) and by the DFG cluster of 
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